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ABSTRACT 

We identify an accretion disk atmosphere and corona from the high resolution X-ray spectrum 
of Hercules X-l, and we determine its detailed physical properties. More than two dozen recom- 
bination emission lines (from Fe XXVI at 1.78 A to N VI at 29.08 A) and Fe Ka,K/3 fluorescence 
lines were detected in a 50 ks observation with the Chandra High-Energy Transmission Grating 
Spectrometer (HETGS). They allow us to measure the density (n e = 2±1 x 10 13 cm from Mg XI), 
temperature {kT = 7 ± 3 eV from Ne IX), spatial distribution, elemental composition, and kine- 
matics (Av < 260 km s -1 ) of the plasma. We exclude HZ Her as the source of the recombination 
emission. We compare accretion disk model atmospheres with the observed spectrum in order 
to constrain the stratification of density and ionization, elemental composition, energetics, and 
thermal stability. The derived disk atmosphere and corona radii are 8 x 10 10 <r<lx 10 11 cm, 
in agreement with previously measured eclipse ingress light curves. The atmospheric spectrum 
observed during the low state is photoionized by the main-on X-ray continuum, indicating that 
the disk is observed edge-on during the low state. We infer the mean number of scatterings (N) 
of Lya and Ly/3 line photons from H-like ions. We derive (N) < 69 for O VIII Lyai, which 
rules out the presence of a mechanism modeled by Sako (2003) to enhance N VII emission via 



-2- 



a line overlap with O VIII. The line optical depth diagnostics are consistent with a flattened 
atmosphere. Our spectral analysis, the disk atmosphere model, and the presence of intense N VII 
and N VI lines (plus N V in the UV) , confirm the over-abundance of nitrogen relative to other 
metals, which was shown to be indicative of CNO cycle processing in a massive progenitor. The 
spectral signatures of a thermal instability in the photoionized plasma are not evident, but the 
measured density is in the stable regime of the models. 

Subject headings: X-rays: binaries — line: formation — line: identification — pulsars: individual 
(Her X-l) — accretion, accretion disks — binaries: eclipsing 

1. Introduction 

High-resolution spectroscopic observations of X-ray binaries have shown in some cases the presence of 
extended X-ray emitting plasmas or outflows surrounding the accretion disk. The detection of these compo- 
nents and their kinematic properties has been made possible by the observation of emission or absorption 
lines from very highly ionized metals. In the case of Hercules X-l, a bright intermediate-mass X-ray binary 
with P or b = 1.7 day orbital period and a P pu ise = 1-24 s X-ray pulsar (Tananbaum et al. 1972; Wilson, 
Scott, & Finger 1997), observations with the XMM-Newton Reflection Grating Spectrometer (RGS) showed 
a strong 35 d phase dependence in the high-resolution spectrum (Jimenez-Garate et al. 2002). The ~ 35 d 
pseudo-periodic cycle (Giacconi et al. 1973) has been associated with a tilted and precessing accretion disk 
(Pcttcrson 1975; Pctterson, Rothschild, & Grubcr 1991; Scott, Leahy, & Wilson 2000). Through the 35 d 
cycle, the X-ray light-curve is asymmetric and contains two maxima: a state of ~ 8 d duration reaching 
the peak flux F max named the main-on, and a secondary high state of ~ 4 d duration reaching ~ 1/3 F max 
named the short-on. A low state with ~ 1/20 F max ensues at other epochs. The main-on spectrum exhibits 
a strong continuum, with evidence for weak and very broad emission lines, while the low state shows bright 
and narrow emission lines dominating over a very weak continuum. The short-on is intermediate, with a 
moderate continuum flux and more intense emission lines than the low state (Jimenez-Garate et al. 2002). 
The Fe Ka line follows a similar pattern. XMM-Newton EPIC data show a 6.4 keV Fe Ka line which is 
practically unresolved during the low and short-on states, and a broad line at 6.5 keV with 330 ± 20 eV 
FWHM during the main-on (Ramsay et al. 2001). Matter close to the magnetosphere and the pulsar appears 
to be observable during the main-on only. This behavior supports the picture of a precessing accretion disk, 
since the spectroscopic appearance of the disk atmosphere and the pulsar should be highly dependent on 
inclination. Analysis of the RGS spectrum showed a clear over-abundance of N with respect to Ne, as well 
as a moderate under-abundance of O and C with respect to Ne. The N over-abundance was proposed as a 
signature of CNO cycle nucleosynthesis occurring at an early phase in the binary system, and it is indicative 
of a star twice or more times more massive than the secondary /companion (Jimenez-Garate et al. 2002). 
This is indicative of a massive progenitor of the system and a period of extreme mass loss, as had been pre- 
dicted by evolutionary models of LMXBs (Podsiadlowski, Rappaport, & Pfahl 2002). Sako (2003) suggested 
that the N over-abundance is not real, but that it is due instead to an X-ray Bowen fluorescence effect which 
occurs at large line optical depths. In this article, we quantify the relevant line optical depths to distinguish 
between the X-ray Bowen effect or the abundance anomaly scenarios. 

Hercules X-l has been observed extensively. Optical light-curves (Gerend & Boynton 1976) and X-ray 
eclipses (Tananbaum et al. 1972) yield a P or b = 1.7 day orbital period. The 1.5 ± 0.3 M Q neutron star has a 
2.3 ± 0.3 Mq companion, HZ Her, which changes from A to B spectral-type over the orbital period, due to 
the strong X-ray illumination on its surface (Reynolds et al. 1997). The unabsorbed luminosity of Her X-l is 
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L = 3.8 x 10 37 ergs" 1 , using a distance of D = 6.6 ± 0.4 kpc (Reynolds et al. 1997). The Her X-l broadband 
X-ray spectrum during the main-on state consists of a blackbody component with temperature kT ~ 90 eV, 
plus a power-law component with a 24 keV exponential cutoff, and a 42 keV cyclotron feature (dal Fiume et 
al. 1998). The X-ray light-curve, the variations in the pulse profiles (Deeter et al. 1991), and the variability 
of the dips (Giacconi et al. 1973; Scott & Leahy 1999), are fit by a geometric model of a processing, warped 
accretion disk, with an 85° inclination with respect to the line of sight, a 20° precession opening angle for 
the outermost disk, and an 11° precession angle for the innermost disk (Scott, Leahy, & Wilson 2000). The 
ultraviolet (UV) spectrum exhibits line emission from C V, N V, and O V, which have two separate velocity 
components. A narrow UV line component is thought to originate on the illuminated face of HZ Her. The 
broad UV line region likely originates in a prograde accretion disk of <~ 10 11 cm radius (Boroson et al. 2000). 
One of the questions to be resolved is whether the X-ray line emission has a similar origin to the UV. The 
broadening of the X-ray line emission indicates that it can originate at radii similar to the UV. For the case 
of an irradiated disk atmosphere, the UV line region is denser and deeper than the X-ray line region. 

In this work, we interpret the broadband (1.5 A to 30 A) high-resolution X-ray spectrum of Her X-l and 
investigate the nature of the line emitting region(s). The Chandra grating spectrum allows us to detect a 
wealth of spectral features which had not been previously observed in Her X-l. Our data reduction procedure 
is described in §2. To interpret the spectrum, we use two complementary approaches: 1) an analytic approach 
for constraining the physical parameters of the plasma, based solely on the photoionized plasma physics; 
and 2) a synthetic approach, by comparing the spectrum with that of an astrophysical model. In §3, discuss 
the origin of the Fe Ka line and we perform spectral diagnostics on density, kinematics, and location from 
the recombination lines. In §4, we utilize standard photoionized plasma models and radiative recombination 
rates to perform an emission measure analysis, and to quantify the ionization distribution of the plasma, the 
emitting volume, and the ratios among the metal abundances. We also compare the model temperatures 
with those measured from radiative recombination continua (RRC). In §5, we quantify the mean number 
of scatterings of selected lines, which shows that the optical depth is not high enough for the X-ray Bowen 
fluorescence effect to take place in Her X-l. Line depth diagnostics also provide evidence for an anisotropic 
plasma distribution. In §6, we identify the nature of the X-ray emission region on the basis of a comparison 
of an accretion disk and corona model with the observed spectrum. In §7, we discuss the thermal stability 
of the plasma. Finally, in §8 we compare the Her X-l spectrum with that of Accretion Disk Corona (ADC) 
sources, and we conclude in §9. 



2. Chandra Observations 

Hercules X-l was observed on 2002 May 5 at 10:15 UT with the Chandra High Energy Transmission 
Grating Spectrometer (Canizares et al. 2000, HETGS). The observation lasted 49.4 ks. Data from the All-Sky 
Monitor (ASM) instrument onboard the Rossi X-ray Timing Explorer (RXTE), indicated that Her X-l was 
in the low state (\& = 0.44 — 0.46) while at mid-orbit (</> = 0.33 — 0.67) during the Chandra observation. From 
the RXTE ASM data, we determined that first main-on start signal immediately preceding the Chandra 
observation occurred on MJD 52384.22. Scott & Leahy (1999) observed that main-on turn-ons occur only 
orbital phases 0.23 or 0.68, possibly due to a disk-orbit locking mechanism. The closest, most likely main-on 
start time occurred at orbital phase 0.23, at MJD 52383.956, ~ 10 hr earlier than the RXTE ASM detection. 

We processed the HETGS data with the CIAO analysis software version 3.0. We obtained Medium 
Energy Grating (MEG) and High Energy Grating (HEG) spectra, which have energy resolution FWHM of 
AA = 0.023 A and AA = 0.012 A, respectively. The ACIS-S CCD was set in FAINT mode since the low 
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state does not produce pile-up in the image, and this minimizes the background in the spectrum. We use 
the destreak tool to clean the ACIS-S image. The mkrmf and mkarf tools are used to produce the response 
matrices and effective area applicable to our observation. We use the calibration data CALDB version 2.26, 
which was updated 2004 Feb 2 to correct for the contamination layer on ACIS-S. We used the ISIS spectral 
analysis software in our model fitting and flux measurements (Houck & Dcnicola 2000). The lightcurve 
routine was used to extract the time variability of the flux. Her X-l was placed at a Y Offset of 0.33 arcmin, 
which is practically on-axis. 



The HEG and MEG spectra of Her X-l exhibit emission lines throughout the X-ray band. Together, 
the spectra exhibit bright lines from the H-like and He-like ions of N, O, Ne, Mg, and Si, plus the lines from 
H-likc ions of S and Fe. The MEG and HEG counts spectra in Figure 1 show a strong Fe Ka fluorescence 
line, as well as numerous recombination emission lines. The HEG spectrum in Figure 2 highlights the Fe 
Ka line from low ionization states and a line from the high ionization state Fe XXVI. The spectra of HEG 
and MEG are added for display purposes in Figure 3, and the MEG spectrum at the high wavelength end 
is shown in Figure 4. The detection of RRCs, as well as the fact that weak Fe L lines are observed, arc 
characteristic of photoionized plasmas (Liedahl & Paerels 1996; Liedahl 1999). We measure a 0.5-7 keV flux 
of Fb.5-7 = 3.3 x 10~ n crgcm~ 2 s _1 , and the fluxes of individual spectral features are listed in Table 1. We 
use Cash statistics to fit the individual spectral features and obtain the statistical errors. Continuum fit 
parameters are shown in Table 2. 



The data allow us to fit independently the Kai and K0J2 lines, which are merged in a single feature. 
We also detect the K(3 line. The measured Fe K fluorescence line energies, as well as the line ratios among 
Kai, Ka2, and K/3, all show that the iron is neutral, based on the Kaastra & Mewe (1993) calculations. 
The Ka/K/3 ratio is 6.2 ± 1.6 (compared to 7.99 for Fe I from Kaastra & Mewe), while the Kai/Ka2 ratio 
is 2.28 ± 0.45 (compared to 2.00, idem). Any iron ion with an L-shcll electron can fluoresce with high yields 
as well. However, the line energy observed corresponds to Fe I-Fe IX (Palmeri et al. 2003), and likely up to 
Fe XIII. We approximate the 6.4 keV fluorescence line flux (in ph cm~ 2 s _1 ) with 



where we take Fe to be the continuum flux (in ph cm~ 2 s _1 keV -1 ) observed during the main-on by dal 
Fiume et al. (1998), rp is the optical depth of M-shell iron, Yp = 0.34 is the fluorescence yield of Fe, f2 is 
the solid angle subtended by the fluorescing plasma from the vantage point of the X-ray pulsar, < T < 1 
is a transmission coefficient due to absorption from species other than Fe I-Fe XIII and Compton scattering, 
and / is the fraction of the fluorescing region which is visible to us. Roughly half the fluorescent photons 
are emitted downward into the optically thick gas. The 0.1 to 100 keV continuum flux during the main-on 
is / Fe<IE = 6.7 x 10~ 9 crgcm~ 2 s _1 , and it is taken from the Model 2 fit by dal Fiume et al. (1998), which 
consists of an absorbed blackbody plus broken power law, with an exponential cutoff and cyclotron absorption 
line. Using tf c ~ 4, which is an estimate of the depth from which an Fe K photon can escape, we obtain 
fTfl ~ 0.22. The solid angle subtended by the companion star is f2 ~ 0.50, and the solid angle subtended by 



3. Spectral Analysis 



3.1. 



Fe fluorescence 
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the disk photosphere is ~ 1.0, as obtained from the illuminated disk model (Jimcnez-Garate, Raymond, & 
Liedahl 2002). These f2 are quite similar. Both cases produce fT ~ 0.2-0.4. Given 0.33 < < 0.67 during 
the observations, we estimate that / ~ 2/3 of the illuminated star is visible. For the disk, f < 1/2 for the 
scenario of a self-shielding edge-on disk (see next section). In both cases reasonable transmission T > 0.5 
can be obtained. Thus, we need to employ the orbital variability to distinguish amongst these components 
(§3.4). The Fe Ka fluorescence line flux and broadening are much larger during the main-on than during 
the short-on, and the low state line flux measured with XMM-Newton, of (2-12) x 10~ 4 ph cm~ 2 s _1 (Zanc 
ct al. 2004), brackets the Fe Ka line flux we measure with HETGS (see Table 1). 

3.2. Unbroadened Lines 

The gratings do not resolve any Doppler width in the emission lines. The Doppler broadening (Av) a is 
measured by fitting Gaussian profiles to the brightest lines in the spectra, such that (Au)fwhm — 2.35(Av) CT . 
Table 1 shows that the best upper limits on Av are placed on the lines with largest A. The observed line 
broadening is entirely attributed to their doublet nature. We set 90 % confidence limits of (AA) ff = 0-8 mA 
for the O VII i line with MEG, (AA) CT = 8-17 mA for the O VIII Lya line with MEG, and (AA) CT = 0-3 mA 
for the Ne IX i line with HEG. The remaining bright lines have similar (AA) CT . The intrinsic line broadening, 
plus any residual instrumental broadening for MEG and HEG not accounted for in the response matrix, 
is quantified using the spectra of a bright star from which no velocity broadening is expected. We choose 
the star HR 1099 as reference, which yields (AA) CT = 4.8 ± 0.4 mA with HEG and (AA) CT = 5.2 ± 0.6 mA 
with MEG at 12.13 A (99% confidence limits). By fitting a plasma model, we verify that this broadening 
can be fully accounted for by the Ne X Lya-doublct, which is separated by 5 mA and should show a 2:1 
intensity ratio. The doublet fully accounts for the measured (AA) ff . Other lines in the HR 1099 spectrum 
are consistent with this. The upper limits on the velocity broadening on Table 1 are based on cither the 
measured or reference (AA) CT , whichever is largest. At the slightly lower resolving power of the XMM- Newton 
RGS, the short-on and low state emission lines of Her X-l were unresolved (Jimenez-Garate et al. 2002). 

The absence of a Doppler velocity broadening constrains the dynamics of the disk atmosphere and the 
geometry of the disk. Since the disk inclination is i ~ 85°, the orbital velocity of the disk is given by 
v ~ yjGMjR cos$, where G is the gravitational constant, M the neutron star mass, R the disk radius, 
and $ is the azimuth on the disk. If the lines produced by a ^-symmetric disk, translating the measured 
(Au)<j into a physical size of the disk requires modeling the line profiles with a calculated (or assumed) 
cmissivity versus radius. The disk lines would have double-peaked profiles. The profiles from a centrally 
illuminated disk of 10 8 5 < r < 10 11 cm radius in full view calculated with the Jimenez-Garate, Raymond, 
& Liedahl (2002) model can be fit with a (Av) a ~ 750 km s _1 Gaussian. Clearly the disk cannot be in 
full view in the Her X-l low state. This is an indication that we are observing the outer rim of the disk 
atmosphere and corona. For a ^-symmetric, edge-on, and optically thin disk, the emission region radius 
would be r > GM/[4(Av) 2 ]. For the X-ray lines observed at highest resolving power, the inferred lower 
limits with this equation are r > (2-3) x 10 11 cm, larger than the neutron star Roche lobe and the disk size 
r ~ 1.4 x 10 11 cm deduced from the eclipse ingress light curves of the He II, C IV, Si IV, and O V emission 
lines in the UV (Chiang 2001; Boroson et al. 2000). Therefore, the disk atmosphere is likely asymmetric. One 
interpretation is that due to the flared disk geometry, the disk shields itself, such that we only observe the 
far side of the disk, which causes the cos$ factor to reduce Av. This does not require an asymmetric disk, 
just one that partially shields itself. Another interpretation is that the disk is asymmetric or warped, which 
can similarly reduce the observed Aw. Evidence of shielding in a warped disk was found by Chiang (2001) 
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in the UV emission line profiles. A third possibility is that the disk is ^-symmetric, but has an atmosphere 
which is orbiting at sub-Keplerian velocities. The case of the asymmetric disk implies a modulation of the 
X-ray line fluxes with orbital phase (see §3.4). 

3.3. A location and density diagnostic from He-like ions 

The Hea line triplets are shown in Fig. 5. At the HETGS resolution, each Hea complex consists of the 
intcrcombination (i), forbidden (/), and resonance (r) lines. The He-like triplet diagnostics were performed 
with the XMM-Newton RGS spectrum for N VI, O VII, and Ne IX (Jimenez-Garate et al. 2002). Both RGS 
and HETGS spectra show that the R = f/i flux ratio is zero for the latter ions. The HETGS spectrum 
reveals that the Mg XI and Si XIII Hea lines behave differently: the R ratio is nonzero and increasing with 
Z (see Table 3). The value of the R ratio depends on the atomic kinetics in each ion. The / line can get 
converted into the i line by a process in which the metastable ls2s 3 Si level is excited to \s2p 3 Pi, 2 (Porquet 
& Dubau 2000). This occurs when the excitation rate Wf^i is larger than the decay rate Wf of ls2s 3 Si to 
the ground state. We note that the particular ls2s 3 Si — > \s2p 3 P2 transition does not modify the R ratio 
for the low and mid-Z ions, for which ls2p 3 P 2 decays back quickly to ls2s 3 Si, instead of decaying to the 
ground state. An R ~ ratio can be produced by either collisional excitation above a critical electron density 
n cnt (p or q ue t & Dubau 2000; Bautista & Kallman 2000) or by photoexcitation above a critical UV-photon 
flux F" lt (Mewe & Schrijver 1978). We discuss both as limiting cases for which R — > 0. If both the plasma 
density and the ambient UV flux are below their thresholds, the R ratio reaches an asymptotic value of 
R > 2.2, which increases with T and depends weakly on Z. For Mg XI at T = 3 x 10 5 K, the maximum 
R = 2.8 (Porquet & Dubau 2000). 

We first discuss the limiting case in which the density is below threshold and the UV flux is above 
threshold. From the measured UV flux F„, we constrain the distance (d) between the X-ray line emission 
region and the UV source. The dominant source of UV flux is the surface of the companion (see below). A 
similar diagnostic was used by Kahn et al. (2001). Jimenez-Garate et al. (2002) set a limit of d < 7 x 10 11 cm 
for Her X-l. With the newly detected lines, d is shown to be equal or larger than the binary separation. To 
measure d, we quantify the photoexcitation rate by using the UV fluxes that were observed with the Hubble 
Space Telescope (HST ) an the Far Ultraviolet Spectroscopic Explorer or FUSE (Boroson et al. 1996, 2000, 
2001), and with the Hopkins Ultraviolet Telescope (Boroson et al. 1997). The relevant photoexcitation rate 
is: 

7Tg 2 

Wf-i = F u if osc , (2) 

m e c 1 

where e, m e are the electron charge and mass, c is the speed of light, / osc is the oscillator strength, and F v , i 
is the flux (in photons s _1 cm~ 2 Hz -1 ) at frequencies resonant with the ls2s 3 Si — > ls2p 3 Pi transition 
in the UV. We use the / osc calculated with the HULLAC atomic code (Klapisch et al. 1977), and those 
calculated by Cann & Thakkar (1992). Since F Vf ^ i is a function of d, we constrain d by relating uif_>; to Wf. 

In the second limiting case, the UV flux is below threshold and the density is above threshold. In this 
case, the lower limit of n e can be obtained from the upper limit on R using the Porquet & Dubau (2000) 
calculations. The d calculated with Eqs. 2 and 3 in the first limiting case, as well as the inferred n e in the 
second limiting case, are shown in Table 3. 

A valuable new constraint is inferred from the Mg XI Hea lines, which show a critical R ratio, that is, 
the value of R = 0.52 ± 0.13 is intermediate between the asymptotic values of R = 0.0 and R = 2.8. The 
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measurcd R ratio implies that 

iOf_>i - tOf . (3) 

For He-like ions other than Mg XI, Eq. 3 becomes an inequality. We interpret the Mg XI R ratio in two 
limits and in the general case. In the low density limit, the UV flux is at a critical value, and F v is measured. 
In the low F v limit, the density is at a critical value, and the measurement of n e is applicable. In the general 
case, the same numbers for Mg XI become a lower limit for d and an upper limit for n e . A confidence region 
in {d,n e ) space requires new atomic kinetics calculations. Since in the general case d > 2 x 10 11 cm, we 
conclude that the X-ray recombination line emission cannot possibly originate in the illuminated atmosphere 
of HZ Her. This does not preclude an Fe Ka fluorescence line to be produced on HZ Her. If the recombination 
lines did originate within the strong UV field of the companion, the Mg XI and Si XIII would have shown no 
/ line and a strong i line, like the rest of the He-like ions. The caveat for the rf-limits is the possibility that 
the UV from HZ Her are shielded on their way to the X-ray emission region, but not on their way to our 
line of sight. The 864 A photons which photoexcite Si XIII are particularly vulnerable since they are above 
the Lyman edge. Such a shielding material could be the accretion disk itself, in which case the origin of the 
recombination lines in the disk atmosphere would be proven as well. Previously, no distinction could be made 
between the case of disk coronal emission and emission from the illuminated companion (Jimenez-Garate et 
al. 2002). 

We base our arguments above in the fact that the primary source of UV emission in the binary system 
is HZ Her and not the accretion disk. This is known from the orbital phase variability. The UV continuum 
increases gradually from <f> ~ 0.2 until it peaks at orbital phase tfi ~ 0.5, and then decreases gradually until 
the illuminated face of HZ Her comes out of view at <p ~ 0.9 (Boroson et al. 2000). The intermediate mass 
(~ 2.3 M Q ) and size of HZ Her causes it to produce copious UV emission by reprocessing of the neutron star 
X-rays. The accretion disk emission in the 1260-1630 A band, which is identifiable from eclipse ingresses and 
egresses at <fi = ±0.1, represents a small fraction (~ 1/10 to 1/20) of the peak HZ Her emission at that band 
(Boroson et al. 2000). However, the question remains of how much UV continuum from the disk is visible 
to the disk atmosphere. Although the observed UV from the disk is fainter than that from the companion, 
that is partly because the disk is seen edgc-on. 

In sum, the weakening of any photoexcitation effects in the R ratios in Mg XI and Si XIII places a lower 
limit on the distance from HZ Her to the X-ray emitting region, or a requirement on shielding from the UV. 
These limits show that the X-ray lines are not produced on the illuminated surface of the companion, but 
instead arise in a high density region that is likely associated with the accretion disk. 

3.4. Search for orbital variability 

The X-ray light curve is flat at mid-orbit. The light curves for the total and Fe Ka line fluxes are 
shown in Figure 6. A Bayesian block analysis (Scarglc 2003) on these light curves showed no deviations from 
constant flux with 99.9% confidence. The coverage is limited to 1/3 of the orbit, so we cannot entirely rule 
out orbital variability in the spectrum. After all, Her X-l exhibits eclipses in the low state (Scott & Leahy 
1999). Clearly, the variability of the Her X-l emission is distinct from that observed in the ADC source 
4U 1822-37, the latter likely arising on an illuminated disk bulge (Cottam et al. 2001). The UV continuum 
emission from HZ Her shows a ~ 25% flux variation from <fi = 0.5 to <j> = 0.65, peaking at mid-orbit, and 
the lines of O V and N VI decrease in strength A(f> = 0.2 away from eclipses (Boroson et al. 2000). With 
better phase coverage than ours, XMM-Newton UV flux and Fe Ka flux measurements show a clear orbital 
modulation (Zane et al. 2004, see their Fig. 9), which suggests that at least some of the Fe Ka flux originates 
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on the face of the companion. Our data show that the Fe Ka flux was unusually flat during our observation. 

4. XSTAR Photoionized Plasma Models and Emission Measure 

The differential emission measure (DEM) maps the ionization distribution of the plasma, ft serves as a 
basis for testing astrophysical emission models and measuring elemental abundance ratios. A DEM analysis 
is valid insofar as the spectrum is dominated by recombination emission, which occurs when the optical 
depth is such that the plasma does not destroy or enhance the lines above their recombination values. The 
luminosity of a line from level u to level I is given by 

dL ul = n e n Zti+1 E ul r] ul a RR dV (4) 

where n z ,i+i is the number density of an ion with atomic number Z and charge +(i + l), E u i is the transition 
energy, dV is the differential volume, a RR = a RR (T) is the total radiative recombination rate in units of 
cm 3 s _1 , and r\ u \ is the fraction of all the recombination photons that produce line emission through the 
transition This can be related to the DEM = J d(E M) / d(\og £_) through the d(EM) = n\dV 

available at each log£, since: 

dL ul = ^±lE ulVul a RR d(EM) = !^±±E ul S ul d(EM) (5) 
n e n e 

where the ionic density n Zj j+i = Aznu fi+i is often expressed in terms of the elemental abundance Az, the 
proton density n Rl and the charge state fraction of the recombining ion /j+i. The specific line power is 
defined as S u i = r] u ia RR . Thus, the total line luminosity is 

T A z n H f d(EM) 

Lui = E u i / fi+iSui-}- — -i-d(log£) (6) 

This sort of DEM procedure was illustrated by Liedahl (1999) and applied by Sako, Licdahl, Kahn, & Paerels 
(1999), and it is analogous to what has been used in collisionally ionized plasmas as d(EM)/dT (Kaastra et 
al. 1996). 

In the simplest analysis, we estimate the total emission measure EM = Jn 2 e dV for each ion, assuming 
each is formed at a single £ Q and T a . The EM for every prominent line in the spectrum is shown in Figure 
7, calculated with the Solar abundances from Grevesse & Sauval (1998), with the updated C and O values 
from Allende Prieto, Lambert, & Asplund (2002). We estimate the elemental abundance ratios by enforcing 
the continuity of the EM as a function of £ Q . The N VII and N VI line fluxes show a N over- abundance with 
respect to C, O, Nc, Mg, Si, S, and Fe. The C VI line flux is obtained from the XMM-Newton RGS spectra 
obtained by Jimenez-Garate et al. (2002). From Equation 6, it is useful to adopt e = fi+\S u i (in units of 
photons cm 3 s _1 ) as a measure of the line emissivity from now on. This e is a function of £. We calculate 
the f i+1 in a grid of models with 0.0 < log£ < 5.0 and A£ = 0.1, using XSTAR (Kallman & McCray 1982). 
We use the main-on spectrum as the ionizing continuum. The S u i are calculated with HULLAC (Klapisch 
et al. 1977). The log T a and log £ of formation for the Lya or Hea lines are defined as the average 
values of log T and log £ weighted by e. The XSTAR/HULLAC model temperatures of formation span 
32, 500 K < T a < 6.7 x 10 6 K, for N VI through Fe XXVI, respectively. We calculate the recombination rate 
at a single £ Q and assume f i+ i = 0.25. We use f i+ i = 0.25 because the log ^-averaged, emissivity- weighted 
value of is approximately equal to that value. While this is just an approximation, it is more realistic 
than simply setting f i+ i = 1. The EM analysis in Figure 7 has the advantages of being simple and of making 
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the result independent of the plasma models; however, one sacrifices the accuracy of the EM distribution 
with £. 

In principle, a more accurate analysis than the above is obtained with DEM = d(EM)/d(log£), by 
assuming that the DEM is constant in the £ range at which each ion is formed. Again, we enforce continuity 
of the DEM as a function of £. We use the aforementioned Solar abundances. Figure 8 shows that the N 
abundance has a large excess, C is moderately depleted, and O is near the Solar value. A recalculation of 
the DEM with the Wilms, Allen, & McCray (2000) Solar abundances establishes that the N excess is just 
as large, while both the C and O depletion appear significant. This stems from the fact that the Wilms et 
al. abundances for both C and O are larger than Allende Prieto et al.'s by ~ 62% and ~ 74%, respectively. 
The abundance pattern, and especially the N abundance, is clearly a signature of CNO cycle products, 
as shown previously by the XMM-Newton RGS spectral analysis, where the DEM was parameterized as a 
power-law and the Wilms ct al. abundances were used (Jimenez-Garate et al. 2002). The uncertainties in the 
Solar abundances do not affect this conclusion, since no significant O depletion is needed. The CNO process 
consists of the rapid conversion of C to N and the 10 3 times slower conversion of O to N (Clayton 1983, 
and references therein). Evolutionary models of 2-12 M Q stars produce an increase in the N abundance by 
factors of a few, with C depletion in the tens of percent and O depletion of a few percent or less, after the 
H-burning phase is over (Schallcr, Schaerer, Meynet, & Maeder 1992). Clearly the most noticeable signature 
will be the N abundance. It can be seen in the DEM that the abundance ratios among Ne, Mg, Si, S, and Fc 
are close to the Solar values. Figure 8 shows that successive pairs of the H-like and He-like ions of a given 
element indicate that the DEM decreases with £, from log£ ~ 1.5 to log£ ~ 2.4. There, the DEM turns over 
and flattens out. The consequences of the DEM shape being rather flat are explored in the next section. 
This DEM is sensitive to the modeled thermal stability of the plasma, as will be discussed in §7. 



4.1. Volume, Scale Height, and Density Diagnostics 

Consider the filling volume V of each line emission region in the optically thin limit. Note that the 
DEM is rather flat, i.e. it varies by less than an order of magnitude (once abundances are accounted for). 
Since £ cx nj 1 , a flat DEM requires that V oc £ 2 . Therefore, the region emitting Fe XXVI fills ~ 10 2 
times more volume than the Si XHI-emitting region, and ~ 10 more volume than the N Vl-emitting region. 
The large disparity in the volume filled by different ions can be interpreted in two ways. The first way is 
consistent with a model in which Fe XXVI is associated with hot corona with ~ 10 4 times the scale height 
of the compact, flattened disk atmosphere which produces N VI. A scale height ratio of the same order was 
calculated from the hydrostatic model atmosphere (Jimenez-Garate, Raymond, & Liedahl 2002). A second 
way is a scenario in which the coldest, densest part of the plasma is concentrated in small clumps, while the 
hottest, more diffuse part of the plasma is distributed in a large volume. 

Combining plasma diagnostics, photoionization balance, and emission measure, we constrain the size 
and scale height of the line-emitting corona, set robust density constraints, and we estimate optical depths. 
For this purpose, we select Fe XXVI, Mg XI, and O VIII. For each of these ions, we plot the density n e as a 
function of distance r from the X-ray source. The emitting plasma is most likely to reside inside the Roche 
lobe of the neutron star. The n e (r) limits drawn for Fe XXVI in Figure 9 are: 

1. An upper limit on the density set by thermal and ionization balance. If the density is too high, the 
plasma would be too cold to emit the line, because a given line can only be emitted at a fixed range of 
£ = L/(n e r 2 ). Thus, n e < L/(^ m i n r 2 ), where £ m j n is defined such that 95% of the line emissivity is at 
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£, > ?min- We take L to be the main-on luminosity, which is reprocessed into line emission during the 
low state. 

2. A lower limit on the density given by the line flux through the emission measure EM and the maximum 
volume U max , such that n e > y/EM/V max = y/3 EM /{A 7rr 3 ). 

Figure 9 includes the density locus for a spherical shell of Thomson depth r ~ 0.05, which is consistent 
with the above limits. Such is the depth required to Compton scatter the main-on continuum and produce 
the continuum observed in the low state. This region of the corona extends to <~ 10 11 cm. The coronal 
size derived from a shallow component of an eclipse ingress is r > 6 x 10 10 cm: the ingress X-ray light 
curves exhibit a steep and fast component and a slow and shallow component, also implying a core of X-ray 
emission concentrated at r < 10 9 cm (Choi et al. 1994). High resolution spectroscopy is allowing us to probe 
the extended region the corona. For the Fe XXVI corona with r ~ 0.05 to be enclosed by the Roche Lobe, 
the scale height must be h/r ~ 0.1, or more conservatively, 0.01 < h/r < 0.3. 

We set stringent density constraints for the Mg XI emission region, taking advantage of the Mg XI R ratio 
being at its critical value (see §3.3 and Fig. 10). We calculate the same n e limits as for Fe XXVI. In addition, 
we calculate the critical n e (in the low UV flux limit), and the critical r found from UV photoexcitation (in 
the low density limit). These are based on the R ratio, as discussed in §3.3. This allows us to compare the 
derived density with the modeled density, for the case where the plasma is inside the Roche lobe. To obtain 
the model density at the Her X-l main-on luminosity we interpolate between two disk atmosphere models 
at L dd and 0.1 L dd, and we find that the density predicted by the atmosphere model (n e = 3 x 10 13 cm~ 3 ) 
is consistent with that derived directly from the R ratio. This provides support to the disk atmosphere 
model, because photoionization balance alone would allow one or two order of magnitude variations on n e , 
depending on the optical depth of the plasma. The agreement in n e also implies that density rather than 
photoexcitation effects are driving the R ratio in Mg XL 

The O VIII diagnostics in Figure 11 allow us to assign an optical depth to a given geometry. This figure 
shows the two n e limits discussed previously. In addition, we plot the density of the plasma for two simple 
geometries: a spherical shell and a flat pill-box or atmosphere. We show two solutions with optical depths 
that are consistent with the density limits. The case of the spherical shell requires larger edge depths (of 
toviii ~ 3) than the case corresponding to a flat atmosphere, for which toviii ~ 0.01 in the face-on direction, 
as would be the case for a warped disk. The values that we will obtain from the optical depth diagnostics 
in §5 will favor the scenario of a flat atmosphere. 

4.2. Abundances from the (N VI + N VII)/(0 VII + O VIII ) line ratio 

Using the photoionized XSTAR and HULLAC models discussed in the beginning of §4, we quantify the 
relationship between the P = (N VI+N VII)/(0 VII+O VIII) line photon flux ratio and the N/O elemental 
abundance ratio. Using the P ratio as an abundance diagnostic has the virtue of being independent of 
the EM and DEM analyses, relying instead on the XSTAR plasma models and the recombination rates 
obtained from the HULLAC atomic code. The XSTAR models are calculated over a grid of log 10 £ values 
that cover the range where the lines of interest are produced. We assume that the line emission is dominated 
by the radiative recombination process, as is shown to be the case for Her X-l in §5. Figure 12 shows 
the calculated E = [e(N VI)+e(N VII)] / [e(0 VII)+e(0 VIII)] line emissivity ratio as a function of £, 
not accounting for elemental abundances (e was defined earlier in §4). Therefore, we have P = E (N/O), 
where (N/O) is the fractional abundance ratio between nitrogen and oxygen. The E and P ratios do not 
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depend on density. At practically any £ at which the lines of interest are produced, 0.7 < E < 2.4. The 
latter values represent the worst case error for one assuming P ~ (N/O). The measured line ratio ranges 
from P = 0.87 ± 0.11 (Jimencz-Garate et al. 2002) to P = 1.54 ± 0.42 (this work). Combining these 
numbers, we obtain (N/O)/(N/O) = 5.9 ± 0.6, using the Allcndc Prieto et al. (2002) and Grevesse & 
Sauval (1998) Solar abundances as reference, and (N/O)/(N/O) = 9.2 ± 1.0 using the Wilms et al. (2000) 
Solar abundances. The error bars on these abundance ratios only include statistical errors. Compare this to 
the 4 < (N/O)/(N/O) < 8 that can be estimated by enforcing continuity on the DEM in Figure 8. 



5. Quantifying optical depth 

X-ray spectroscopic analysis allows us to measure or set limits on the mean number of scatters of 
line photons, or to the line optical depth. The photoelectric optical depth can also be constrained by the 
spectrum, despite the absence of absorption edges. These measurements can have important implications for 
the geometry of the plasma. They also allow us to validate the previously determined elemental abundance 
ratios. 



5.1. The mean number of scatterings for Lya and Ly/3 photons 

The mean number of scatterings is constrained by the observed Lya/Ly/3 ratios. The Lya/Ly/3 ratio is 
a temperature diagnostic in collisional plasmas. In contrast, this ratio depends very weakly on temperature 
for a recombining plasma, as we have verified from the HULLAC recombination rates. A Ly/3 photon from 
a H-like ion can get resonantly absorbed by another ion of the same species (see Fig. 13). There is a 
probability 7 = 0.12 that the absorbed Ly/3 photon is converted to photons with lower energy (including 
an Ha). Otherwise, the absorbed Ly/3 is re-emitted. The probability of a Ly/3 photon to survive N line 
scatterings is: 

P = (1 - 7 ) W = 0.88 w , (7) 

and therefore the Ly/3 photons cannot scatter more than a few times before getting destroyed. As the 
number of scatterings increases, the Lya/Ly/3 intensity ratio increases with it, because Lya photons do not 
get destroyed in line scatterings (7 = 0). None of the observed Lya/Ly/3 line ratios shown in Table 4 show 
evidence for multiple line scatterings, except for O VIII. Since we deal with a large ensemble of photons, the 
number related to the line ratio is the mean number of scatterings (N), instead of N above. 

The observed Lya/Ly/3 ratio allows us to constrain the mean number of scatterings of Ly/3 photons 
directly, and of Lya photons by inference. If 7 > 0, the mean number of scatterings (N) for a slab geometry 
satisfies (Hummer & Kunasz 1980) : 

w = ^t < 8 » 

where El is the line energy lost to resonant absorption and Eq is the total line energy generated in the gas 
column. There are two relevant cases for the upper limit of (N). In the case 7 = 0, photons do not get 
destroyed (as for Lya), then (N) does not have an upper bound, and Equation 8 does not apply. In the case 
7 > (as for Ly/3), (N) is bounded above (by {N) < 7.33). as E L /E G -> 1 in Equation 8. 

We first use Equation 8 for O VIII. We detected 7 counts of the O VIII Ly/3 line in the 15.995 A 
to 16.015 A range. Given the continuum and background level of 3.0 counts, this represents a detection 
significance of 95%, as derived from a Poisson distribution. The caveat is that this is only valid if the 
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background is uniform in A, and it may be subject to systematics. Since this is not a line search (i.e. the 
line energy is fixed), this may be a positive detection, but it needs confirmation. Taking Ly/3 at face value, 
we derive E L /E G < 0.9 for O VIII Ly/3 with 90% confidence, and obtain (N) < 6.6. The detected line is 
centroided at 16.005 A (O VIII Ly/3), and it is one MEG-resolution FWHM away from Fe XVIII at 16.023 A. 

To infer the mean number of scatterings of Lyai from those of Ly/3, we need to relate (N) to r, the 
mean line optical depth. Hummer & Kunasz (1980) performed such a calculation for a slab geometry with 
a uniform source of line photons, 7 = 0, and Voigt parameters of 4.7 x 10~ 2 to 4.7 x 10~ 4 . They obtained 
(N)/t ~ 0.55-0.6 for r = 1, of (N)/t ~ 0.60-0.65 for t = 10, and (N)/t ~ 0.73-0.85 for r = 100. The 
Voigt parameter is the ratio between the natural broadening and the velocity broadening of the line. Define 
c T = (N)/t. To generalize this result for 7 > 0, note from Equation 2.15 in Hummer & Kunasz (1980) that 
(1 — 7) can be factored out of the expression for (N), such that (N) = c t t(1 — 7) is applicable in that case, 
with c T and r as before. Note c T = c T (r) is a function of r. 

The mean line optical depth t is proportional to the oscillator strength r oc / osc . Therefore, we obtain 
(N) oc c T / osc (l — 7), from which wc estimate (N) <~ 30 for O VIII Lyai photons. A similar number can be 
obtained through the more simplistic Equation 7. We ignored the continuum optical depth. The effect of a 
nonzero continuum depth is to decrease (N) below the derived values. Up to this point, this result depends 
on the detection of the Ly/3 line, but there is another way to deduce (N) for O VIII Lycui photons. 

We set an upper limit on (N) for O VIII Lyai without relying on O VIII Ly/3, by using the Ne X lines. 
This is valid since both ions are produced at a largely overlapping ionization parameter range. We find that 
(N) < 2 for Ne X Ly/3 photons (at 90% confidence). Using the c T once again, this implies (N) < 9 for 
Ne X Lyai. The EM derived from Ne X and O VIII are the within 12% of each other, consistent with being 
equal given the statistical uncertainties. Scaling with the Wilms, Allen, & McCray (2000) or Allcnde Prieto, 
Lambert, & Asplund (2002) Solar abundances, we deduce (N) < 69 for O VIII Lyai. Scaling with CNO 
abundances, (N) is smaller than or equal to the latter limit. Both are consistent with the (N) ~ 30 value 
for O VIII Lyai which was derived above. 

In contrast, the Lya/Ly/3 ratios of Mg XII and Si XIV on Table 4 may show evidence for enhancement of 
Ly/3 due to photoexcitation effects at moderate column densities. Photoexcitation of ions by X-ray continuum 
photons, followed by radiative decay, can dominate the emission mechanism. This process, also referred to as 
resonant scattering, can produce a characteristic Lya/Ly/3 ratio (Kinkhabwala et al. 2002). At low column 
densities, both lines are unsaturated, resonant scattering dominates over recombination emission, and the 
Lya/Ly/3 ratio simply depends on the oscillator strengths. At moderate column density, the Lya/Ly/3 ratio 
is reduced below the low column value, when the Lya line reaches saturation but the Ly/3 line does not. The 
evidence for photoexcitation of Mg XII and Si XIV at the moderate column density regime is a ~ 2cr result, 
so we do not discuss this scenario further. At high column densities, recombination emission dominates over 
resonant scattering, because resonant scattering has completely saturated while the emission measure keeps 
growing with the column, and the Lya/Ly/3 ratio is determined by the ratio of the recombination efficiency 
of the lines, which happens to be similar to the oscillator strength ratio. In this section we have expanded 
upon the Kinkhabwala et al. (2002) picture, showing that at high line optical depths, the Lya/Ly/3 ratio 
can increase above its pure recombination value due to the destruction of the Ly/3 line. 
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5.2. The optical depth of the r line in He-like ions 

The value of the G = (i + f)/r ratio depends on whether the plasma is collisionally ionized or pho- 
toionized, but it can also can be used to quantify the mean optical depth of the r line. The G ratio depends 
on optical depth for three reasons: 1) the r line has much larger oscillator strength than either the i or / 
lines; 2) the recombination rates of the i, / lines are larger than those of the r line, because the i, f lines 
have larger statistical weight than the r line; and 3) the G ratio depends very weakly on temperature in a 
photoionized plasma. This diagnostic has been quantified by Kinkhabwala et al. (2002) and Wojdowski et 
al. (2003). In the limit of small line optical depth, we have G ~ 0.3, because the r line is enhanced relative 
to the / and i lines by continuum X-rays that photo-excite the ground state (Wojdowski et al. 2003). In 
the small line depth regime, the emission is dominated by resonant scattering over recombination emission. 
In this regime, the line spectrum is dependent on the spectral energy distribution. In the limit of large line 
optical depth, G reaches an asymptotic value of 4.5 (Porquet & Dubau 2000), which depends very weakly 
on the the T values derived from photoionization equilibrium, and recombination emission dominates over 
resonant scattering. 

The observed G ratios shown in Table 3 are consistent with the limit of large line optical depth. By using 
the lower bounds on G, we can set lower limits on the mean depth of the r lines as defined by Wojdowski 
et al. (2003), of r r > 100 for G > 4 (the maximum value applicable to Ne IX and O VII), r r > 20 for G > 3 
(applicable to Mg XI), and r r > 10 for G > 2 (applicable to Si XIII), assuming a continuum power-law index 
of a = 1.0. The observed main-on index is a = 0.85, but G depends weakly on a. 

5.3. An X-ray atmosphere illuminated at a small grazing angle has anisotropic line depths 

If the plasma geometry is anisotropic, the line optical depth will depend on direction. A first special 
direction vector corresponds to the path of the X-ray continuum photons from the X-ray pulsar, since this 
radiation is photoionizing the plasma. The second special direction is that which minimizes the optical 
depth, since photons will tend to escape from the plasma that way. Based on this, we define two relevant 
depths: an "illumination" optical depth Tj, and an "escape" optical depth r e , which are well-defined and not 
equal if the plasma distribution is illuminated along a direction where the optical depth is not minimized. 
The illumination depth n determines how much of the continuum radiation is absorbed, and therefore how 
much is reprocessed into resonant and recombination emission lines. On the other hand, once a line photon 
is produced, it is more likely to escape in the direction where the optical depth is minimized than through 
other directions, and this corresponds to the escape depth r e . 

The Hea G ratio is sensitive to n, because G depends on the ratio between the resonant scattering and 
recombination line fluxes, as mentioned in the previous section. Further scattering of the r line after it is 
produced does not destroy it unless the continuum depth is significant. On the other hand, the Lya/Ly/? 
ratio is sensitive to both Tj and the escape depth r e . Again, resonant line absorption of Lya or Ly/3 produces 
a dependency on n, because both photons resonant scatter strongly. In addition, the Ly/3 photon can 
be destroyed and converted to Lya if it scatters too many times inside the plasma, and this is why the 
dependence on the escape depth r e arises. 

The flat, geometrically thin shape of the model disk atmosphere is anisotropic, and the atmosphere is 
illuminated by the neutron star at a very shallow grazing angle that was determined from self-consistent 
calculations (Jimenez-Garate, Raymond, & Liedahl 2002). This geometry is shown schematically in Figure 
14. In the atmosphere model, the calculated grazing angle is such that illumination depth is 25 times 
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larger than the escape depth at r = 10 11 cm, and therefore the observed Hea r line depth is expected to 
be much larger than the Ly/3 depth. Sample vertical column densities in the disk atmosphere model are 
Novm = 5 x 10 17 cuT 2 , JVovn = 2 x 10 17 cm~ 2 , and JV Si xiii = 1-6 x 10 16 cm" 2 . These columns imply 
escape line depths of r e ~ 46 (150 km s~ 1 /v) for O VIII Lyai, r e ~ 53 (150 km s~ 1 /v) for O VII r, and of 
r e ~ 1.4 (150 km s _1 /w) for Si XIII r, where v a fiducial velocity width that includes thermal and turbulent 
components. Setting v to the thermal value at kT = 25 eV for O VIII yields r e = 178 for Lyai, and setting 
kT = 3.5 eV for O VII yields r e = 544 for He r. Note that tj ~ 25t c for the disk model geometry, i.e. for 
Si XIII we get n = 36 (150 km s~ 1 /v). These are the expected line depths from the disk atmosphere model, 
modulo the assumptions made about the gas kinematics. 

Compare the above values with the optical depth measurements from Lya/Ly/3 that yielded (N) ~ 30 
for O VIII Lyai, which is equivalent to t ~ 46 given the calculated (N)/t ratios and Voigt parameters from 
§5.1. This optical depth matches the model at the chosen 150kms _1 velocity width. If the O VIII Ly/3 line 
detection is not real, from Ne X it follows that for O VIII Lyai still r < 100. This would simply require a 
smaller velocity width than above, so there is parameter space available to achieve consistency of the data 
with the disk atmosphere model in this case. For the case of O VII, an upper limit r r > 100 was obtained 
from the G ratio, which is easily accommodated by r t ~ 1300 (150 km s~ 1 /v) for the illumination depth. 
This large leaves plenty of parameter space to be consistent with the other lower limits on r r set for 
Ne IX, Mg XI, and Si XIII. For example, Si XIII was measured to have r r > 10, which is consistent with the 
atmosphere model tj = 36. For reference, the model escape depth for the continuum is r e (OVIII) = 0.05, an 
illumination depth Ti(OVIII) = 1.2. 

Our results above depend on the estimated kinematics, and therefore are only approximate consistency 
checks. Since the velocity field in the disk is produced by both Keplerian and turbulent motions, r e will in 
fact be highly anisotropic, and it will depend on the photon momentum vector relative to the disk velocity 
field. To self-consistently derive a turbulent velocity from the Lya/Ly/3 ratio, one requires a calculation of 
the line transfer which accounts for this velocity field. 

5.4. RRC/RR ratios and electron temperature 

The electron temperature measured from the Ne IX RRC is kT = 7 ± 3 eV, or T = 81000 ± 35000 K. 
This was measured by fitting the semi- Maxwellian profile of the Ne IX RRC at A = 10.37 A, the brightest 
RRC in the spectrum, detected at the 5a level (see Fig. 3). The XSTAR plasma code and the HULLAC 
recombination rates yield emissivity- weighted average temperatures kT a which are shown on Tables 4 and 5. 
The T Q for Ne IX is in good agreement with the measured T. The significance of detection for other RRCs 
obtained with the Cash statistic ranges from the 2a to the 5a level (see Table 1). To fit the RRCs, we rebin 
the spectra to bin sizes A A of 0.05 A to 0.09 A (coarser than the bin size used for fitting lines). While most 
RRCs do not have sufficient statistics for an accurate T measurement, T is left as a free parameter and the 
best-fit values are consistent with T Q for most cases. 

In Table 5, we show the observed RRC to RR (radiative recombination lines) flux ratios, compared to 
the HULLAC calculations, in the recombination-dominated case, for both a single temperature model and 
for the disk atmosphere model. The RRC/RR ratios for O VII, O VIII, and Ne X are significantly weaker 
than expected for the single-T model, but agreement is generally better for the disk atmosphere model, which 
accounts for a broad distribution of temperatures 2 < kT < 860 eV. It is therefore likely that the RRCs are 
composed of multiple temperatures, with the low-T component being observed more readily than the high-T 
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component. The high-T RRC component, which usually dominates the recombination flux, is too broad to 
be separable from the continuum. In this favored scenario, the weak RRCs are simply an indication of the 
multi-T nature of the plasma. 

5.5. Verification of CNO abundances: the line optical depths are not large enough to 

convert O VIII Lya to N VII Lya 

We test for the presence of the X-ray Bowen fluorescence mechanism by which a O VIII Lya line becomes 
a N VII Lya line, as calculated by Sako (2003). We determine whether this X-ray Bowen fluorescence effect 
can account self-consistently for the apparent excess intensity observed in N VII Lya, N VI i and r, and the 
N VII RRC with respect to other lines, such as the O VIII Lya and O VII i and r. We find that: 

1. The X-ray Bowen fluorescence mechanism is unlikely to enhance N VI Hea emission in the same 
proportion as N VII Lya. This would require line opacities to overlap in energy space under nearly 
identical conditions as those found for N VII Ly( and O VIII Lya2, whose energies are 1 Doppler 
width away at T — 50 eV. We checked for coincidences between the energies of the brightest lines of 
abundant elements in the 22.46 to 26 A band and the resonance line energies of N VI. We did not 
find line overlaps closer than 46 Doppler widths at T = 50 eV, and the number of Doppler widths is 
likely to be higher since models show a peak N VI Hea flux at T = 2.8 eV. As shown in the HETGS 
spectrum in Fig. 4, as well as from the XMM-Newton RGS spectrum (Jimenez-Garate et al. 2002), 
both N VI Hea and N VII Lya indicate a N over-abundance. 

2. The measured mean number of scatterings of O VIII Lya and Ly/5 photons is two orders of magnitude 
lower than required for the X-ray Bowen fluorescence effect to take place. The measured line ratio 
between O VIII Lya and Ly/3 depends on (N), due to the destruction of Ly/3 or its conversion to Lya. 
We derive (N) ~ 30 for O VIII Lyai. This was verified via the same ratio in Nc X, which scales to 
(N) < 69 for O VIII Lyai with Solar or CNO abundances (see §5.1). Compare this to the (N) ~ 5 x 10 3 
scatterings needed for the N VII enhancement (Sako 2003). 

3. The detection of the N VII RRC with XMM-Newton and Chandra provides additional evidence for 
a N/O overabundance. At the escape line depths implied by the measured Lya/Ly/3 ratios, the 
RRC fluxes are unaffected by optical depth effects. The N VII RRC detected with XMM-Newton RGS 
(Jimenez-Garate et al. 2002) at 10cr has flux (6.5±2.8) x 10~ 5 ph cm~ 2 s _1 , with 90% confidence errors, 
and HETGS detected it at 3.5er with flux (18 ± 5) x 10~ 5 ph cm^s" 1 , with 68% confidence errors 
(Table 1). These fluxes are larger than the < 1 x 10~ 5 ph cm -2 s -1 expected from Solar abundances 
and the fluxes of O VII RRC and O VIII RRC. 

4. An analysis of the UV line emission from Her X-l using the same sort of illuminated disk model 
concluded that N was over-abundant in HZ Her, with a N/C abundance ratio of ~ 2 (Raymond 1993), 
which implies a N overabundance relative to Solar of ~ 6. This shows that both UV and X-ray disk 
emission models yield consistent results. 

5. The velocity shear in the disk atmosphere and corona increases the line escape probability, which 
would decrease (N) below the static values calculated by Sako (2003) for a given column density. This 
suppresses X-ray Bowen Fluorescence. 
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We conclude that the large nitrogen over-abundance interpretation of the Her X-l spectrum by Jimenez- 
Garate et al. (2002) is valid and cannot be self-consistently explained by the X-ray Bowen fluorescence effect. 
The large nitrogen over-abundance implies that significant CNO processing occurred in a massive star in the 
system. New transfer calculations are required to improve the optical depth and dynamical constraints on 
the accretion flow. Such calculations need to include thermal and ionization balance, stratification, and gas 
orbital kinematics in order to calculate line transfer more accurately than existing models. 

6. Comparison with an accretion disk atmosphere and corona model 

In Figs. 2-4 we over-plot the synthetic spectrum of an illuminated accretion disk atmosphere and corona 
model from Jimenez-Garate, Raymond, & Liedahl (2002). We compare the data to this disk atmosphere 
model. The model assumes that the neutron star continuum energizes the vertically stratified structure of a 
hydrostatic accretion disk in thermal and ionization balance. We used two adjustable parameters to fit the 
spectrum: a normalization factor of N — 0.11, and the disk radius range 10 10 9 < r < 10 11 cm (the outer 
annulus in the model). Both parameters affect the normalization and do not influence the spectral shape. 
We did not include any Doppler velocity broadening effects. The ionizing continuum in the model consists 
of a bremsstrahlung with T — 8 keV and a luminosity L x = 10 38 3 (D/6.6 kpc) 2 TV ergs -1 . Given the 
approximate linear dependence of line flux on continuum flux (Jimenez-Garate, Raymond, & Liedahl 2002), 
the effective ionizing luminosity in the model fit was L x — 2.3 x 10 37 ergs -1 . This is nearly the unabsorbed 
luminosity during main-on state of Her X-l in the 0.1-200 keV band, L = 3.8 x 10 37 ergs -1 (dal Fiume et 
al. 1998). The model normalization can be interpreted as simply a geometrical parameter which indicates 
that the observed emission is produced by an outer section in the accretion disk atmosphere and corona with 
area A <~ 3 x 10 21 cm 2 . That A is a small fraction of the total disk area is in agreement with observing an 
outer disk rim, as suggested by the small velocity broadening of the emission lines (§3.2). The neutron star 
continuum observed during the main-on powers the recombination emission observed during the low state. 
The precession of the accretion disk allows us to isolate the reprocessed emission from the direct emission 
from the neutron star. 

The relative strengths of the spectral lines, which are fixed in the model, are a direct measure of the 
structure of the accretion disk atmosphere and corona. In general, they constrain the temperature, ionization, 
and density distributions. In the context of our model, the relative line strengths are indicators of the 
vertical ionization structure of the disk atmosphere and corona, since the ionization level and temperature 
are stratified with vertical depth. In Table 1, we compare the relative strengths of lines with different 
ionization parameters (except for Fe Ka fluorescence at 1.94 A, which is fit separately, because it is not yet 
calculated in the disk model). We note that the relative strengths of the Lya lines of S XVI, Si XIV, Mg XII, 
Ne X, and O VIII lines were fit with ~ 50% or better accuracy, with virtually zero degrees of freedom, because 
N and the radius range, the two free parameters in the model, only affect the overall flux normalization. 
Since O VIII and S XVI are produced at ionization parameters that are one order of magnitude apart, it is 
a notable success for the model to be able to reproduce their flux ratios. In the disk model context, this 
indicates that the emission measure distribution of the upper half of the atmosphere agrees well with the data. 
The incidence angle on the disk atmosphere is self-consistently determined by the model (Jimenez-Garate, 
Raymond, & Liedahl 2002). 

The Fe XXVI Lya line probes the hottest and most extended region of the disk corona observed with 
HETGS (see §4.1 and Fig. 2). The XSTAR photoionization equilibrium plasma model (Kallman & McCray 
1982), coupled with the recombination rates obtained with the HULL AC atomic code (Klapisch et al. 1977), 
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show that the electron temperature of the plasma producing the Fe XXVI line is T e ~ 6.7 x 10 6 K, with 
ionization parameter log£ ~ 3.8. Similar parameters are obtained with the disk atmosphere model. The 
observed Fe XXVI Lya line flux is 3.9 times larger than the disk model's. Judging from the disk model, the 
Fe XXVI Ha and H/3 recombination lines at 9.68 A and 9.53 A (see Fig. 3 and Table 1) are 10 times weaker 
than expected. The flux ratio of Ha and Lya lines is not expected to depend strongly on temperature in 
the pure photoionization case. The Fe XXVI Lya line was also detected with XMM-Newton (Zane et al. 
2004). An emission line of Fe XXV at 1.85 A appears in the model, as well as a Fe XXIV complex with lines 
at 11.18 A, 11.03 A, and 10.62 A (th ese lines appear in Fig. 3). None of these relatively weak lines are 
detected. In short, the upper coronal structure appears to be more ionized and to have a larger emission 
measure in Her X-l than in the disk atmosphere and corona model. 

We have shown that several key characteristics of the X-ray recombination spectrum can be explained 
by an outer section of the accretion disk which is being illuminated by the neutron star. This disk at- 
mosphere model describes much of the Her X-l X-ray line spectrum from first principles. However, there 
are discrepancies between the model and the data. One notable discrepancy is N VII and N VI, which are 
under-predicted by a factor of ~ 10. This is consistent with the §4.2 results. Also, the He-like ion line fluxes 
are systematically underpredicted by the model, as we discuss below. 

7. Thermal stability of the photoionized plasma and the H-like/He-like ion line ratios 

A thermal instability has been predicted to be present in plasmas which arc photoionized by an X-ray 
continuum. The "type I" instability is produced by an imbalance between bremsstrahlung cooling and Comp- 
ton heating, just below the Compton temperature (Krolik, McKee, & Tarter 1981). The "type II" instability 
is produced by the imbalance between recombination cooling and photoionization heating (Kallman 1984). 
The photoionized disk atmosphere models predict a type II instability regime at 6 x 10 4 < T < 5 x 10 5 K 
(Jimenez-Garate, Raymond, & Liedahl 2002). X-ray recombination emission is suppressed in the latter 
regime. The recombination line fluxes are sensitive to the ion fraction The instability reduces the 

interval [log£ Q — A(log£), log£ + A(log£)] at which f i+1 > 0.1. For this reason, the disk atmosphere model 
predicts very weak Mg XI line fluxes relative to Mg XII, contrary to observations (as shown in Fig. 3). The 
DEM analysis indicates that the Mg XI flux is above the expectation (Fig. 8). The ions with £ Q near the 
instability regime will have their line power reduced and their derived DEM increased. 

The He-like ion lines of Si XIII, Mg XI, Ne IX, O VII, and N VI are underpredicted by the atmosphere 
model. Also, the DEMs derived from the same He-like ions are systematically larger than the DEMs from 
the H-like ions. This may imply that the plasma is more thermally stable than the models predict, or it 
may require a larger EM of T ~ 10 5 K plasma. The He-like ion lines are especially sensitive to the thermal 
instability because of a key physical difference: the H-like ions form by the recombination of fully stripped 
atoms, so their line emissivities are less sensitive to £ than those of He-like ions. The for the He- like 

lines are single-peaked functions, while the for the H-like lines are asymptotic curves with f i+1 — ► 1 

as £ — + oo, and f i+1 — > as £ — > 0. This causes the line emissivities for the He-like ions to have narrow 
peaks and small A(log£), while in contrast, the line emissivities of H-like ions have broad asymmetric peaks 
with large A(log£), with tails that extend to large £. This is why, in spite of a possible overlap of the £ of 
formation between the He-like and H-like ions, the H-like ion line fluxes will still dominate over the He-like 
ion lines at high £. Therefore we expect the He-like ion lines to be more sensitive indicators of the ionization 
distribution and the instability regime than the H-like ion lines. The changing trends of the EM in Figure 7 
compared to the DEM in Figure 8, suggest that the thermal instability regime is larger in the models than 
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in the plasma, but this is difficult to ascertain without a systematic study of the plasma models. 

The 7i e measurement from Mg XI Hea (§4.1) also tests stability of the photoionized plasma, since the 
disk atmosphere model predicts a regime of forbidden density. The n e regions allowed by both the spectral 
diagnostics and disk atmosphere models are shown in Figure 10. The density of formation for Mg XI is 
n e = 3 x 10 13 cm~ 3 according to the disk atmosphere model. The forbidden densities are just below that, 
at 2 x 10 12 < n e < 3 x 10 13 cm -3 . The measured n e = (2 ± 1) x 10 13 cm~ 3 is located in the high density 
end of the stable regime. The thermal instability calculation in the disk model is in agreement with the 
observations in this case. 

The X-ray spectroscopic data is now sufficiently detailed that it is possible to start testing the regimes of 
thermal instability in the photoionized plasma. On the one hand, a DEM analysis based on XSTAR plasma 
models and the HULLAC recombination rates suggests that the the plasma is more thermally stable than 
was originally predicted. On the other hand, the density constraints from the R ratio of Si XIII indicate that 
the density is just at the edge of the stability region, as derived from the disk atmosphere model. Further, 
the same disk atmosphere model predicts the Mg XI Hea flux to be suppressed due to the instability, but 
no such suppression is observed. 

8. X-ray Spectral Comparison with other LMXBs 

We compare the spectral properties of the Her X-l low state with two low-mass X-ray binaries: 4U 1822- 
37 and 2S 0921-63. Both 4U 1822-37 and 2S 0921-63 are classified as ADC sources. These LMXBs exhibit 
X-ray emission lines with unresolved velocity broadening. The limits on the velocity widths (FWHM) of 
the recombination lines measured with Chandra HETGS are Av < 1500 km s _1 for 2S 0921-63 (Kallman, 
Angelini, Boroson, & Cottam 2003), Av < 300kms _1 for 4U 1822-37 (Cottam et al. 2001), and Av < 
260kms~ 1 for Her X-l (from O VII in Table 1). Her X-l and 4U 1822-37 exhibit RRCs that indicate 
that photoionization is the dominant ionizing mechanism in the plasma. In 2S 0921-63, the bright O VII 
intercombination line (giving the value of G ~ 4 for that diagnostic line ratio) is indicative of a photoionized 
plasma (Kallman, Angelini, Boroson, & Cottam 2003). The non-detection of RRCs in 2S 0921-63 is likely 
due to the low signal-to-noise ratio of the emission features. The weak or undetected Fe L-shell lines in 
these LMXB spectra are yet another signature of photoionization. It is therefore clear that the lines in these 
LMXBs arc produced by recombination following photoionization. 

When looked at more carefully, a comparison of the spectra of 4U 1822-37, 2S 0921-63, and Her X-l 
reveals a nuanced picture. First, the continuum level varies significantly from source to source: 4U 1822-37 
has ~ 10 times larger continuum flux in the 0.5-7 keV band than 2S 0921-63 or Her X-l. The continuum 
in 4U 1822-37 is harder than in the other two sources. Second, the soft X-ray line spectra in 4U 1822-37 
and Her X-l appear to be quite similar, i.e. the fluxes of the Si XIV, Mg XII, Mg XI, Ne X, Ne IX, O VIII, 
and O VII lines are roughly within tens of percent of each other, albeit the exact line ratios and density 
diagnostics differ. On the other hand, 2S 0921-63 shows weaker line fluxes than 4U 1822-37 and Her X-l, 
perhaps because the eclipse is blocking the source of the line emission. Third, the line ratios between Fe 
emission lines from both high and low ionization species are different. Her X-l exhibits a much more intense 
Fe Ka line from Fe I-Fe XIII at 1.94 A than 2S 0921-63, but a nearly equal flux of the Fe XXVI line at 1.78 A. 
2S 0921-63 exhibits both Fe XXV and Fe XXVI lines with nearly equal fluxes, and a Fe Ka fluorescence line 
with 1/4 of the flux of the Fe XXV/Fe XXVI lines (Kallman, Angelini, Boroson, & Cottam 2003). It is 
conceivable that the fluorescing material was occulted by the eclipse of 2S 0921-63, while the hotter, more 
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extended emission was not. 4U 1822-37 also showed Fe Ka and Fe XXVI lines, with Fe Ka being brighter 
(Cottam et al. 2001), but it did not show the large flux contrast between those lines observed in Her X-l. 
Fourth, the N VII Lya and N VI Hea lines are most intense in Her X-l (see Fig. 4). So far, other LMXBs 
and ADC have not shown such strong N lines, albeit the measurement is rendered difficult for most cases 
due to the large absorbing hydrogen columns Ah that are common in these systems. In §4 and §5.5, we 
quantified this abundance anomaly, which was also observed with the XMM-Newton RGS (Jimenez-Garate 
et al. 2002). 

We now compare X-ray flux variability with orbital phase. While our observations of Her X-l are 
restricted to mid-orbit at 0.33 < 4> < 0.67, the 2S 0921-63 observations were restricted to mid-eclipse, and 
the 4U 1822-37 observations spanned full orbits. The emission during the 2S 0921-63 eclipse originates in 
either the outer disk rim or in an extended circumsource region (such as a corona) with an upper size limit 
of 7 x 10 10 cm (Kallman, Angelini, Boroson, & Cottam 2003). The emission line region in 4U 1822-37 
was identified to be the illuminated disk bulge, because the line fluxes were observed to peak at 4> ~ 0.25 
(Cottam et al. 2001). Flux variations were not observed in Her X-l, but with the limited phase coverage we 
can ascertain that it does not show the same sort of variation than 4U 1822-37. Judging from the spectral 
variations as a function of orbital phase, the azimuthal location of the line emission regions in Her X-l, 
2S 0921-63, and 4U 1822-37 appears to be distinct amongst these sources. 

Her X-l, 4U 1822-37 and 2S 0921-63 are sources that are observed at high-inclination angles (85°< 
i < 90°), which causes the accretion disk to block the direct X-ray emission from the accreting object. 
ADCs have a higher optical-to-X-ray flux ratio than other LMXBs, which is a well-known indicator that the 
accreting object is not observed directly (White, Nagase, & Parmar 1995, and references therein). This, in 
turn, causes the EWs of the X-ray emission lines in ADCs to be much larger than in other LMXBs, because 
only ~ 10 -2 -10 -3 of the available continuum energy is reprocessed into recombination lines. This scenario 
is consistent with the model presented in §6, where only a fraction of the accretion disk area can account 
for the line emission. Incidentally, the only other ADC for which a high-resolution spectrum was obtained 
is AC 211 (in the M15 cluster), but it did not show emission lines (White & Angelini 2001), probably due 
to contamination from a bright and nearby LMXB burster. 

In sum, high-inclination LMXBs exhibit recombination-dominated line spectra and unresolved line 
broadening. The illuminated disk bulge seen in 4U 1822-37, however, is not observed in Her X-l with 
the current data, and only the extended emission during eclipse was observed in 2S 0921-63. A definitive 
comparison of the accretion flow geometry of these sources will require new spectral data with good orbital 
phase coverage for Her X-l and 2S 0921-63. 

9. Conclusions 

1. We detect more than two dozen spectral features in the X-ray spectrum of the Hercules X-l low state 
with the Chandra HETGS. Most of the observed features correspond to the recombination signatures 
expected from a centrally illuminated accretion disk atmosphere and corona located in the outer ac- 
cretion disk, given the following evidence: 

• Using photoexcitation diagnostics with the He-like ion lines, coupled with the observed ultraviolet 
continuum fluxes by Boroson et al. (2001), we ruled out the possibility that the recombination 
emission lines originate on the illuminated companion star, HZ Her. 

• The spectral model for a centrally illuminated accretion disk atmosphere and corona fit most of 
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the observed fluxes and ratios of the X-ray recombination lines. The model requires the low state 
emission to be energized by the continuum observed in the high state. The caveat is that the 
He-like lines are underpredicted by the model, and the observed line broadening is smaller than 
expected for a plasma corotating with the accretion disk. The small Doppler widths may be due 
to vignetting effects that allow only the far edge of the disk to be visible, as would occur for an 
edge-on disk with a flared or warped geometry. In support of this, the model fits indicate that 
only the outer rim of the disk (8 x 10 10 < r < 10 11 cm) is observed, in contrast with LMXBs 
such as EXO 0748-67, where much of the inner disk emits observable X-ray lines ( Jimenez-Garate, 
Schulz, & Marshall 2003). 

• The measured R line ratio of Mg XI allows us to precisely determine a density of n e = (2 ± 1) x 
10 13 cm~ 3 , which agrees with the disk atmosphere model predictions (n e = 3x 10 13 cm~ 3 ). This 
corresponds to the low UV field limit, which is the most likely scenario for Mg XI. Using UV 
observations (Boroson et al. 2001), we deduce that the disk may be far enough from HZ Her for 
the Mg XI triplet to be unaffected by photoexcitation. The presence of UV emission and high-n e 
necessitates the re-calculation of R ratios in (n e , i*uv) space. 

• The X-ray spectrum probes radically different length scales in the atmosphere and corona. For 
example, an emission measure analysis shows that the volume where Fe XXVI is emitted must 
be ~ 10 4 times larger than that for N VI. Compare this to the cool model atmosphere emitting 
N VII at T ~ 32500 K, with 10 7 cm thickness, and 3 x 10 13 cm -3 density, which coexists with a 
10 10 cm thick hot corona emitting Fe XXVI, with 10 n -10 12 cm" 3 density and T ~ 10 6 K (with 
intermediate layers of ionization in between). 

• The RRC/RR ratios and the measured T = 81000 ± 35000 K from Ne IX RRC are generally in 
agreement with the disk atmosphere model. Single-T models do not reproduce well the observed 
RRC/RR flux ratios of H-likc ions, while the disk model produces improved agreement. This 
indicates the presence of multiple temperature components, since any high-T component of the 
RRCs is lost in the continuum. 

• The line depths derived from the spectrum via the Lya/Ly/3 and G ratios are consistent with the 
ionic column densities and flattened geometry of the model disk atmosphere, where the first ratio 
is most sensitive to escape depth (vertical depth) and the second is determined by the grazing 
incidence depth. To obtain consistency, we use a turbulent velocity width of Av ~ 150 km s -1 . 
We rely on the relationship between r and the mean number of scatters (N) calculated for a range 
of Voigt parameters by Hummer & Kunasz (1980). 

2. We find a prominent Fe K«i,a 2 , K/3 complex indicative of a neutral plasma. We cannot ascertain 
whether the Fe K fluorescence is produced in the disk or on HZ Her, or both. The Fe Ka line flux 
from HETGS does not vary with orbital phase at 0.3 < </> < 0.7, but Zane et al. (2004) did observe a 
clear modulation with <f> with XMM-Newton. 

3. The lines in the X-ray spectrum have unresolved velocity broadening (Av < 260 km s -1 ), similar to 
what has been observed for the lines in ADC sources. This is evidence that the accretion disk is edge-on 
with respect to our line of sight. Previously observed spectral changes with the XMM-Newton RGS 
showed a transition between broad lines in the main-on state and narrow lines in the short-on and low 
states, which suggested emission from the inner disk during main-on and from the outer disk during 
short-on and low state (Jimenez-Garate et al. 2002). The well-known 35 d flux modulation has long 
been attributed to the precession of an edge-on disk. 
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4. By measuring the mean number of line scatterings, we rule out the presence of the X-ray Bowen 
fluorescence effect calculated by Sako (2003). We constrain the mean number of scatterings of the 
O VIII Lyai line to (N) < 69. The upper limit is based on the detection of O VIII and Ne X Ly/3 lines, 
which must be destroyed after a large number of scatterings. In contrast, the X-ray Bowen fluorescence 
effect would require (N) ~5x 10 3 . Therefore, the N VII/O VIII Lya line ratio is unaffected by optical 
depth effects. 

5. The signature of CNO abundances observed with XMM-Newton (Jimenez-Garate et al. 2002) is con- 
firmed. Strong N VI Hea, N VII Lya, and N VII RRC lines are present in the HETGS and in the XMM- 
Newton RGS spectra. The N abundance is verified to be higher than Solar. The simplest EM analysis, 
the DEM analysis, and the (N VII + N VI)/(0 VII + O VIII) line ratio yield 4 < (N/O)/(N/O) < 9. 
The uncertainty is dominated by systematics in the analysis rather than by the effective area calibra- 
tion. The depletion of carbon derived from the RGS data is also evidenced in our DEM analysis. This 
is supported by measurements of the N/C ratio made by comparing the observed UV spectrum to that 
from an illuminated accretion disk model (Raymond 1993). The relative O abundance is sensitive to 
the Solar abundance values, being under-abundant with the Wilms, Allen, & McCray (2000) values, 
and consistent with Solar with the Allende Prieto, Lambert, & Asplund (2002) values. 

6. Given the problems with fitting the He-like ion line fluxes with cither the disk model or the DEM, we 
deem that either a) the thermal instability regime is over-predicted; or 2) the disk atmosphere model 
requires more emission measure at T ~ 10 4 -10 5 K. The He-like ion lines are more sensitive probes of the 
thermal equilibrium of the plasma than the H-like ion lines. It is possible that the thermal instability 
has been over-estimated in the models, therefore cutting the He-like ion fluxes disproportionately. 
However, the density diagnostics seem to be consistent with the calculated stable density values. 

7. There are line transfer effects which remain to be accounted for in the models, such as the competition 
between line emission by photoexcitation and line destruction, and the effects of repeated photoioniza- 
tions and recombinations produced by lines and RRC. However, the Her X-l spectrum indicates that 
these are second-order effects, compared to the pure recombination emission. 

8. The atmospheric and coronal components unveiled by the X-ray line spectra are more extended (10 10 — 
10 11 cm), cooler, less luminous, and have smaller vertical Thomson depth (tt ~ 0.01) than the compact 
(10 9 cm) ADC component with tt > 1, which resides close to the compact object. Indeed, both compact 
and extended ADC components are evident in the eclipse ingress light curves observed in Hercules X-l 
(Choi ct al. 1994). Therefore the plasma observed in high resolution spectra is not part of what is 
commonly referred to as an ADC (i.e. a fully ionized spherical distribution of plasma that is Thomson 
scattering the neutron star X-rays), but instead may be a continuation of it at the Roche lobe scale. 
The extended ADC component observed in eclipse coincides with our size measurements for the X-ray 
line region. 
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Fig. 1. — The MEG and HEG count spectra of Hercules X-l during its low-state shows numerous radiative 
recombination lines and continua, a strong Fe Ka fluorescence line, and a continuum which may be produced 
by Compton scattering. The large range of temperatures at which these emission lines originate is a direct 
probe of the ionization structure of the upper layers of the accretion disk atmosphere and corona. 
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Fig. 2. — HEG spectrum, showing both Fe I-Fe XIII Ka and K/3 fluorescence lines, as well as the Fe XXVI 
Lya radiative recombination line. The disk corona and atmosphere model is over-plotted in red, except for 
the fluorescence lines, which were fit independently. The green label indicates that the line was produced in 
the model but not detected. 
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Fig. 3. — Added HEG and MEG spectra. The disk corona and atmosphere model is over-plotted in red. 
The green label indicates that the line was produced in the model but not detected. 
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Fig. 4. — MEG spectra which show that the nitrogen lines are quite bright relative to the oxygen lines, as 
observed with XMM-Newton (Jimenez-Garate et al. 2002). The disk corona and atmosphere model is shown 
in red. 
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Fig. 5.— He-like triplet lines for O VII, Ne IX, Mg XI, and Si XIII. The line ratios evolve smoothly from O VII 
and Ne IX with R ~ 0, to Si XIII with R = 3 ± 1. The R ratio of Mg XI has a critical value R = 0.51 ± 0.12, 
which can provide a precise measure of the density. The Gaussian line fits (not from the disk model) are 
over-plotted in red. 
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Fig. 6.— Light curves of: a) the total HETGS dispersed count rate with A > 1.5 A, and b) the Fe I-Fe XIII 
Ka line HETGS count rate. The orbital period is 1.7 days. 
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Fig. 7. — Estimated emission measure (EM) at a single ionization parameter of formation £ OJ derived from 
the HETGS line fluxes (diamonds) and from the C VI flux measured with RGS (Jimenez-Garate et al. 2002, 
filled circle). Both N VII and N VI indicate an excess nitrogen abundance. This rough estimate of the EM 
has the virtue of being insensitive to thermal instabilities. We assume a charge state fraction f i+ i — 0.25 
for the recombining ions. 
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Fig. 8. — Differential emission measure (DEM) vs. ionization parameter. The line fluxes measured with 
HETGS (diamonds) are complemented with the fluxes from XMM-Newton RGS (Jimenez-Garate et al. 2002, 
filled circles). This is more accurate than Fig. 7 because the DEM is calculated in a grid spanning the full 
£ range. However, it is also sensitive to the thermal and ionization balance solution from XSTAR (Kallman 
& McCray 1982), as well as the DEM slope. The over-abundance of N and depletion of C is observed. The 
predicted thermal instability regime occurs at 1.9 < log 10 £ < 2.4. The excess DEM for the He-like ions in 
this regime suggests that the plasma is more stable than predicted. 
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Fig. 9. — Density limits derived for the Fe XXVI emission region, versus its maximum distance (r) to the 
neutron star. The Thomson depth is plotted for a spherical shell with r m j n = r/3 and r max = r. 
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Fig. 10. — Limits on the density of the Mg XI emission region, versus its maximum distance to the neutron 
star. The R ratio strongly limits the plasma to be close to the density and photoexcitation boundary. The 
irradiated disk atmosphere model prediction is shown by an X, and the thermally stable density regimes in 
that model are shaded in yellow. 
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Fig. 11. — Limits on the density of the O VIII emission region, versus its maximum distance to the neutron 
star. The velocity broadening suggests that the emission region should be close to the Roche Lobe of the 
neutron star. The O VIII edge optical depths for two geometries are shown: for a filled spherical shell with 
?~min = r/3 and r max = r ("shell"), and for a face-on flattened disk geometry ("flat"). 
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Fig. 12. — Line emissivity ratio (E) built from the N VI, N VII, O VII, and O VIII line emissivities obtained 
from the XSTAR and HULLAC models, as a function of ionization parameter log 10 £ . The observed line 
ratio will differ from this emissivity ratio by a factor of An/Ao, the ratio of elemental abundances. At the 
top of the figure, we show the log 10 £ range at which the line emissivities are > 20% of the peak emissivity, 
with the log 10 £ Q of formation (peak emissivity) marked. 




Fig. 13. — Abridged energy level diagram for a H-like ion, showing the decay channels of the 3p 2 P?,/2,i/2 
levels. These channels show how a Ly/3 photon gets scattered or converted to other photons. Once a Ly/3 
photon is absorbed, there is an 88% probability for an electron at either of the 3p 2 Pz/2,i/2 levels to re-emit 
a Ly/3 photon, resulting in a resonant scattering event. Alternatively, there is a 12% probability for the 
same initial states to decay by emission of an Ha photon. Once an Ha is emitted, the metastable 2s 2 5 , 1 / 2 
level decays to the ground state via two photons (not shown), unless the electron density is above a critical 
threshold n cr jt- For n e > n cr it, the 2s 2 S\/2 level gets excited to 2p 2 Pz/2, and a Lya photon is emitted. The 
photons shown in the figure are all unresolved doublets. 
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Fig. 14. — In an accretion disk atmosphere which is illuminated with 9 <C 1, two optical depths are relevant. 
In the Her X-l spectrum, recombination emission dominates over resonance scattering, implying that > 100 
for the He- like ion r lines. The mean number of scatterings and any line destruction are sensitive to r e . The 
X-ray Bowen fluorescence effect is excluded because (N) < 69 and thus r e < 99 for O VIII Lyai. The 
measured t are beginning to show that the emission region is elongated. 
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Table 1. X-ray Emission Features and Disk Atmosphere Model Predictions 



Linc(s) 


A [A] a (At)) 


a [km/s] 


ni J TIM b A /T „ ,] 

Ubscrvcd r lux Moo 


el Flux b (A) 


N VI i 


29.08 




64 ± 22 


1.9 


N VII Lya 


24.78 




25 ± 4.5 


2.0 


O VII / 


22.10 




1.2 ± 1.2 (< 3.9) 


0.36 


O VII i 


21.800 ± 0.006 


< 110 


40.5 ± 5.5 


9.2 


O VII r 


21.60 




3.1 ± 1.6 


2.5 


O VIII Lya 


18.970 ± 0.006 


< 270 


13 ± 2 


9.0 


N VII RRC 


18.59 




18 ± 5 


0.63 


Fo XVII 


17.05, 17.10 




< 3.8 


0.56 


O VII RRC 


16.78 




2.5 ± 1.0 


4.2 


O VIII Ly/3 


16.01 




0.85 ± 0.30 


2.0 


O VIII RRC 


14.23 




5.8 ± 1.2 


2.4 


No IX / 


13.70 




0.17 ± 0.17 (< 0.96) 


1.4 


No IX i 


13.553 ± 0.003 


< 160 


7.3 ± 0.75 


0.71 


No IX r 


13.45 




0.87 ± 0.33 


0.66 


No X Lya 


12.135 ± 0.003 


< 240 


3.2 ± 0.4 


1.8 


No IX RRC 


10.37 




1.82 ± 0.36 


0.80 


No X Ly/3 


10.24 




0.59 ± 0.16 


0.67 


Fc XXVI Ha 


9.53, 9.67 




< 0.24 


0.79 


Mg XI / 


9.31 




0.84 ± 0.17 


0.32 


Mg XI i 


9.229 ± 0.003 


< 380 


1.65 ± 0.21 


0.14 


Mg XI r 


9.17 




0.59 ± 0.15 


0.10 


No X RRC 


9.10 




0.73 ± 0.19 


0.32 


Mg XII Lya 


8.417 ± 0.003 




1.12 ± 0.28 


1.1 


Mg XII Ly/3 


7.11 




0.42 ± 0.11 


0.30 


Si XIII / 


6.74 




1.08 ± 0.15 


.38 


Si XIII i 


6.69 




0.36 ± 0.12 


.18 


Si XIII r 


6.647 ± 0.003 




0.54 ± 0.13 


.11 


Si XIV Lya 


6.185 ± 0.003 


< 340 


1.44 ± 0.23 


1.8 


Si XIV Ly/3 


5.22 




0.61 ± 0.26 


0.35 


S XVI Lya 


4.73 




1.19 ± 0.31 


0.88 


Fo II-XIII Ka 2 


1.940 ± 0.003° 


< 680 


10.61 ± 1.84 




Fc II-XIII Kai 


1.936 ± 0.003° 


< 680 


22.00 ± 1.94 




Fc XXVI Lya 


1.78 




7.37 ± 1.30 


1.9 


Fc II-XIII K/3 


1.76 




5.10 ± 1.20 





Note. — The model corresponds to an accretion disk atmosphere and corona. The statistical 
errors and upper limits are given to 66 % confidence, except for the upper limits, which are 90 % 
confidence. Symbols: A — line wavelength; (Av) CT — standard deviation of the Gaussian emission 
line. We include systematics in the measurements, except for the line fluxes. We compare the data 
to the disk atmosphere and corona model. 

a Wcak features are assigned nominal wavelengths. Otherwise, the measured wavelengths arc 
shown with error bars. The nominal wavelenghts are rounded up to 0.01 mA. 

b In units of 10 — 5 photons cm" 2 s _1 , 

c This is the relative shift of the Fc Ka pair with wavelength and flux ratios tied. 
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Table 2: Continuum fit parameters 



Component (Model) 


Parameter (Unit) 


Value 


Absorption (1) 


N H (10 18 cm" 2 ) 


13 (fixed) 


Power- law (1) 


Norm, at 1 keV (phot keV -1 cm~ 2 s -1 ) 


(9.6 ±0.4) x 10~ 4 


Power- law (1) 


Photon Index 


0.27 ±0.03 


Blackbody (1) 


Norm. (10 39 ergs _1 / (10 kpc) 2 ) 


(3.5 ±0.2) x 10~ 5 


Blackbody (1) 


Temperature (keV) 


0.181 ±0.015 


(1) 


x 2 /dof 


1292/1844 = 0.70 


Absorption (2) 


Nh (10 1S cm" 2 ) 


13 (fixed) 


Power- law (2) 


Norm, at 1 keV (phot keV -1 cm~ 2 s -1 ) 


(1.34 ±0.01) x 10~ 3 


Power-law (2) 


Photon Index 


0.53 ± 0.01 


(2) 


x 2 /dof 


1703/1844 = 0.92 



Note. — DOF = degrees of freedom. 



Table 3: Helium-like line diagnostics and photoexcitation of the ls2s 3 Si level 





R = f/i 


G = (f + i)/r 




Flux F Xf _ ri 


tfjf 


W{^i(r— 10 li cm) 


d 


n e 


Ion 


Line 


Line 




(10~ 13 erg cm -2 


2 3 Si -> 1 


2 3 Si -. 2 3 Po,i, 2 


Radius 


Density 




Ratio 


Ratio 




s- 1 A- 1 ) 


Rate (s- 1 )' 2 ' 5 ' 


Rate (s _1 ) 


(cm) 


(cm- 3 ) 


Si XIII 


3.0 ± 1.0 


2.7 ± 0.7 


864 


~ i.i^> 


3.56 x 10 5 


7 X 10 4 


> 8 X 10 1U 


< 5 X 10 12 


Mg XI 


0.51 ± 0.12 


4.2 ± 1.2 


1033 


~3.1< 4 > 


7.24 x 10 4 


4 x 10 5 


2 X 10 11 


(2 ± 1) x 10 13 


No IX 


< 0.13 


8.6 ± 3.4 


1270 


1.9 ± 0.1 (3) 


1.09 x 10 4 


5 x 10 s 


< 6 x 10 11 


> 1 x 10 13 


O VII 


< 0.10 


13 ± 7 


1637 


1.5 ± 0.2 (3) 


1.04 x 10 3 


1 x 10 6 


< 3 X 10 12 


> 8 x 10 11 



Note. — The statistical errors are calculated to 66% confidence and the limits to 90%. Symbols: Af = wavelengths of 

ls2s 3 Si — > ls2p 3 Po,l,2 transitions; Wf = radiative decay rates for ls2p 3 Si — > Is 2 1 So ; tff^i = calculated photoexcitation 
rate for ls2p 3 5i — > ls2p 3 Po,l,2- The energy level notation is abridged in the table header. 



References. — (1) Porquet et al. 2001; (2) Drake 1971; (3) from the Hubble GHRS at <f> = 0.56-0.60 by Boroson et al. 1996; 
(4) Boroson et al. 2001, and private communication 2003, FUSE data near <j> = 0.5; (5) From HULLAC atomic code, Klapisch 
et al. (1977). 
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Table 4. Lya to Ly/3 Line Ratios 



Ion 


kT 


Lya/Ly/3 


Lya/Ly/3 




[eV] 


theory 


observed 


VIII 


disk 


4.7 


15 ±6 




25 


4.9 




NeX 


disk 


3.3 


5.4 ± 1.6 




90 


4.6 




Mg XII 


disk 


4.8 


2.7 ± 1.0 




154 


4.8 




Si XIV 


disk 


5.5 


2.4 ± 1.1 




237 


4.8 





Note. — The theoretical values were 
calculated with HULLAC in the optically 
thin limit at the temperature ol peak 
emission, and for the disk atmosphere 
(with temperatures 2 < kT < 860 eV), 
which includes possible line blends with 
other species. T Q ~ average temperature 
weighted by the emissivity of the RR line 
(fi+l S u i). 
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Tabic 5. RRC to Radiative Recombination (RR) Line Ratios 



Ion 




RRC/RR 


RRC/RR 




[eV] 


theory 


observed 


N VI 


2.8 


0.24 






disk 


0.26 




N VII 


14 


0.83 


0.72 ±0.25, 0.28±0.14W 




disk 


0.32 




VII 


3.5 


0.32 


0.06 ±0.02 




disk 


0.35 




VIII 


25 


0.87 


0.45 ±0.12 




disk 


0.27 




Ne IX 


6.5 


0.25 


0.22 ±0.05 




disk 


0.29 




NcX 


90 


1.02 


0.23 ±0.07 




disk 


0.18 





Note. — The theoretical values were calculated with 
HULL AC in the optically thin limit. For RR, we use ei- 
ther Lya or Hea, corresponding to each H-likc or He-like 
ion, respectively. T = average temperature weighted by 
the emissivity of the RR line. The disk model tempera- 
ture range is 2 < kT < 860 eV. 

References. — (1) Jimenez-Garate et al. (2001). 



